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A gaseous Bose-Einstein condensate (BEC) irradiated by a far off-resonance laser has long-range 
interatomic correlations caused by laser-induced dipole-dipolc interactions. These correlations, 
which are tunable via the laser intensity and frequency, can produce a 'roton' minimum in the 
excitation spectrum — behavior reminiscent of the strongly correlated superfluid liquid helium II. 
PACS: 03.75.Kk, 32.80.Qk, 34.20. Cf, 67.40.Db 



According to the celebrated Bijl-Feynman for- 
mula £Q for the excitation spectrum of helium II 



the peculiar "roton" minimum at k s=s 2ir/ro, 
where tq is the average atomic separation, is due 
to a corresponding peak in the static structure 
factor S(k) = {0\p k p< k \0}/N. Here N is the num- 
ber of atoms of mass m, |0) the ground state of 
the system, and f>k = J2 q Cq^q+k the density fluc- 
tuation operator. S(k) is the Fourier transform 
of the pair correlation function and hence pro- 
vides a measure of the degree of pair (2nd order) 
correlation between the atoms. The existence of 
strong pair correlations in helium II may at first 
seem surprising since it remains a liquid even 
at temperatures approaching absolute zero pre- 
cisely because of weak interatomic interactions 
(in combination with a small atomic mass) 
However, despite their apparent weakness these 



interactions are very effective because the den- 
sity of the liquid state is such that the average 
atomic separation, tq = 4.44A, is close to the 
minimum of the attractive interatomic potential 
well at 3 A. 

Contrast this now with an ultra-cold alkali 
atom gas in which the Bose-Einstein condensed 
fraction can be very nearly 100 % 0] . The inter- 
actions in Rb, for example, are repulsive and 
characterized by an s-wave scattering length, 
a w 5.5nm. This is between one and two orders 
of magnitude smaller than the average atomic 
spacing at typical densities. Steinhauer et al [5] 
recently measured the bulk excitation spectrum 
of a 87 Rb BEC and found excellent agreement 
with Bogoliubov theory |3| (appropriate for a de- 
generate almost ideal Bose gas). There was no 
roton minimum, a consequence of the dilutcness 
with respect to a. Indeed, since Eq. (JIJ becomes 
an equality within the Bogoliubov theory £Q one 
sees the pair correlation is small compared to he- 
lium II. Significant pair correlation might exist 
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in gaseous BECs at the very small scale of a, 
but this is fairly inaccessible in such a delicate 
system. 

A marvellous feature of atoms though, is that 
their interactions can be manipulated using ex- 
ternal fields, allowing us to microscopically engi- 
neer the macroscopic properties of a many-body 
system. Thus the experiment of Inouye et al 
[S] took advantage of a Feshbach resonance to 
change the s-wave scattering length using mag- 
netic fields. We have recently proposed the 
use of off-resonant lasers to induce long-range 
dipole-dipole interactions whose characteristic 
length is the laser wavelength. These interac- 
tions can cause laser induced self- "gravity" in a 
BEC, leading to 3-dimensional self-trapping and 
electrostriction accompanied by unusual excita- 
tion spectra [7j, as well as "supersolid" struc- 
tures The aim here is to explore how the 
excitation spectrum and, by virtue of (Q, the 
correlations of a gaseous BEC are changed when 
the interatomic potential is modified via laser- 
induced dipole-dipole interactions. This task 
requires a knowledge of the Fourier transform 
(FT) of the total interatomic potential, so we 
turn to this first. 

Consider a BEC confined by a potential 
#tra P = ^oj^x 2 +y 2 ) + ^j^z 2 into a very elon- 
gated cigar shape (oj r » uj z ), irradiated by a 
far off-resonance plane-wave laser (Fig. ^ in- 
set). The laser polarization is along the long 
z-axis of the condensate to suppress collective 
( "superradiant" ) Rayleigh scattering or co- 
herent atomic recoil lasing that are forbid- 
den along the direction of polarization. The far 
off-resonance condition, together with the small 
extent of the BEC along the laser propagation 
direction, enables us to treat the electromagnetic 
field inside the BEC in the Born approximation 
(field at each point is the sum of the incident plus 
once-scattered fields). Then the dipole-dipole 
potential between two atoms of separation r, in- 
duced by far-off resonance light of intensity /, 
wave- vector kL = &Ly (along the y-axis), and 
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Figure 1: The total (s-wave+dipole-dipole) 
ID interatomic potential U^ ot {z) (FT of ©), 
for w r — 3.5Al. A repulsive contact term 
4£ , R a(fc L w r )~ 2 (l + 4I/3)5(ki J z) is not shown. 
Inset: The laser beam and condensate geome- 
try. 

polarization e = z (along the z-axis) is jllj 



U dd (r) 



la 2 (uj) kj 

47TC£n 



Vzz (kL,r) cos ,(k h y). (2) 



Here a(uj) is the isotropic, dynamic, polariz- 
ability of the atoms at frequency uj — ck^ = 
27tc/Al- The pre- factor can be expressed in 
terms of the Rayleigh scattering rate, 7r, as 
Ia 2 kl/(ATrcel) = (3/2)^. V zz is the com- 
ponent of the retarded dipole-dipole interaction 
tensor generated by the linearly z-polarized laser 
light 



V z . 
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fc3r 3 



(l — 3 cos 2 9) ( cos ki,r + k\jr sin k^r) 

- sin 2 9 k\r 2 cosk L r (3) 



9 being the angle between the interatomic axis 
and the z-axis. The far-zone (k^r 3> 1) be- 
havior of J5J) along the z-axis is proportional to 
— sin(fcLr)/(fcL?') 2 and many atoms (400 at den- 
sities of 8 x 10 14 atoms/cm 3 ) may lie within the 
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characteristic interaction volume (A 3 ) of this at- 
tractive long-range potential. As for the elec- 
tron gas and charged Bose gas, mean-field (here 
Bogoliubov) theory applies in this high density 
regime [T2] . 

The laser (dynamically) induced dipole-dipole 
potential is distinguished from the static field 
(r~ 3 ) case ^H^] by a longer range and a huge 
enhancement of atomic polarizability around a 
resonance. For example, in 5 87 Rb atoms are 
magnetically trapped in the maximally stretched 
\5s 2 S 1/2 ,F = 2, M = 2} state. A laser po- 
larized along z is then 7r-polarized and only 
AM = dipole transitions are allowed. If 
the light is detuned by, say, 5 = 27rx(6.5GHz) 
(i.e. 1134 natural line widths) below the Dl line 
(795. Onm) then only virtual transitions to the 
\5p 2 P\/2,F = 2,M = 2) state need be consid- 
ered. We calculate a « 5.0 x l(T 35 Cm 2 /V ( c f. 
the static value 5.3 x l(T 39 Cm 2 /V). 

In terms of the condensate density n(r) at 
zero temperature, we account for atom-atom 
interactions using a mean-field energy func- 
tional of the form H dd + H s , where H dd — 
(1/2) J rc(r)Z7 dd (r-r')n(r')cZ 3 rdV, an d H s = 

(l/2)(4irah 2 /m) J n(r) 2 d 3 r is due to short- 
range interactions, which are described, as is 
usual, by a delta function pseudo-potential (we 
take here the repulsive case for which a > 0). 
By working with the bare dipole-dipole inter- 
action we assume the Born approximation also 
for atom-atom scattering by this long-range part 
of the total potential. We note that the short- 
range (static) part of the laser-induced dipole- 
dipole interaction can cause a shift in a. For the 
laser intensities and detunings considered here 
this shift is small according to existing estimates 

m 

In a radially tight trap it is reasonable to as- 
sume a cylindrically symmetric ansatz for the 
density profile of radial width w r : n(r) = N 
(irw 2 ) -1 n z (z)exp [—(a; 2 + y 2 )/w 2 ], where N is 
the total number of atoms and n z (z) is normal- 
ized to 1 and left general. Denoting the FT 
of the atomic density by n(k) = J exp[— ik • 



r]n(r)d 3 r, then we have H dd = (l/2)(2vr)" 3 
/ U dd (k)h(k)n(~k)d 3 k, where the FT of the 
dipole-dipole potential J2J, U dd (k) = J exp[— ik- 
r] U dd (r) d 3 r, is the real part of 



U dd (k) 
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2 e 2 c \ k 2 x + {k y - k h ) 2 + k 2 
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h 2 - h 2 
n z /t L 



fc2 + (ky + fc L ) 



VI] 



The principal value of the radial integration 
in H dd can be evaluated analytically so that 
the dipole-dipole energy reduces to a one di- 
mensional functional along the axial direction 
H dd = (N 2 /2) $ n z {z)n z {z')U^ d {z - z')dzdz' 

^{N 2 /4n) f^ z (k z )^(-k z )uf d (k z )dk z , where 
n z (k z ) is the FT of the axial density n z (z). The 
one-dimensional (ID) axial potential that ap- 
pears in this expression has the form 



Q(w r ,k z ) = - 
" (fc L uv) 2j 



- — — Q(w r ,k z ) , 

47T£g C 



3=0 



» { E J+1 



(k 2 -k 2 )w 2 



where $l{Ej [z] } is the real part of the generalized 
exponential integral [TS] . The FT of the total (s- 
wave plus dipole-dipole) ID reduced interatomic 
potential is 

Ufotikz) =4E R a({k L w r )- 2 +lQ{w r ,k z )) (6) 

where En — h 2 k 2 i /2m is the photon recoil en- 
ergy of an atom and X is the dimensionless 'in- 
tensity' parameter 



Ia 2 (w)m 
87reHc?i 2 a 



(7) 



It is emphasized that the radial degree of free- 
dom is contained in (jSJ) via the radius w r . The 
coordinate space potential, U z ot (z), is shown in 
Fig. [U 



(4) 



(5) 
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We now have the essential ingredients to com- 
pute the excitation spectrum of the BEC as it 
is the FT of the effective interatomic interaction 
potential that appears in the Bogoliubov disper- 
sion formula Since the influence of radial 
excitations upon the low-energy spectrum can 
be largely frozen out under tight radial confine- 
ment, we shall consider only axial phononic ex- 
citations and assume that the system is infinite 
along this z direction. In terms of the phonon 
momentum p z = hk z , the axial Bogoliubov spec- 
trum is (cf. Eq. Q) 

E B = y/ clpl + (p 2 /2m) 2 = p 2 J (2mS (k z )) (8) 

where c 2 — Trn(0)w 2 U{; ot (k z )/m. n(0) is the 
central density in the cigar, so 7m(0)u> 2 is the 
linear density along the cigar. For the linear 
parts of the spectrum, c z can be interpreted 
as the speed of sound in the gas. Shining 



n(Q) = 8 x I0 20 atoms/m 3 and radius w r = 
3.5Al = 2.78/im, a 'roton' minimum appears 
when X > 0.051 (i.e. I > 0.506W/cm 2 ), al- 
though the dispersion relation is considerably 
altered far before this. The change in the dis- 
persion relation could be observed using Bragg 
spectroscopy as performed in [SJ- Fig. [2] plots 
the Bogoliubov dispersion for X = 0.057 (/ = 
0.565W/cm 2 ). 

Local to the 'roton' minimum at k = k TO t on 
one can write E = A + h 2 (k — fc rot0 n) 2 /2m*, 
and for the parameters above with I = 0.057 
one finds m* = 0.06m. He II has m* = 0.16m 
0. The static structure factor is plotted in 
Fig. El The peak in S(k z ) corresponds to the 
minimum in the energy spectrum. The model 
described here predicts that when X > 0.066 
(/ > 0.654W/cm 2 ) the minimum touches the 
zero energy axis. At this point the system is 
unstable to a periodic, supersolid-like, density 
modulation piTH]. 




Figure 2: The Bogoliubov dispersion relation for 
87 Rb with w r = 3.5A L and n(0) = 8 x 10 20 
atoms/m 3 . For pure s-wave scattering (X = 0) 
the inverse healing length l/£o = y/8Tran(0) = 
1.32fc L . 

a 795. Onm laser upon a 87 Rb BEC of density 



Figure 3: The static structure factor, S(k z ) for 
various laser intensities. Same paramaters as 
Fig. □ 

The laser induced dipole-dipole potential can 
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lead to electrostriction (compression) of a con- 
densate .7 . In the present regime of low laser 
intensity/large detuning the electrostriction is 
negligible (on a scale set by the collapse thresh- 
old J = 3/2 [7JEI)- This regime also ensures 
the absence of two-body bound states in the ID 
reduced potential shown in Fig. ^ a necessary 
condition for the validity of the Born approxima- 
tion for atom-atom scattering by this potential. 
Only when X > 1.3 do bound states appear. 

The interaction © arises from the forward 
scattering of laser photons by atom pairs. At 
large detunings there are two main compet- 
ing processes that can heat a dense gas : A) 
Light-induced transfer of pairs of colliding 
atoms to a quasi-molecular excited state 
followed by dissociation, releasing w US into the 
kinetic energy of the pair [TJ|. This is a den- 
sity dependent effect whose rate can therefore 
be high. Even when the laser is red-detuned 
from an atomic resonance, when two atoms col- 
lide the energy separation between the ground 
state and a molecular excited state (—C^/r 3 ) 
comes into resonance at small distances. How- 
ever, by choosing 5 so that the resonance point 
occurs between two molecular vibrational states 
this process is suppressed [TJ\. Below the Dl 
line there are only discrete molecular vibrational 
states (i.e. no continuum states) so a detuning 
can be selected which is between these molecular 
resonances |18j . which are narrow at ultra-cold 
temperatures. 

B) Incoherent light scattering by single 
atoms occurs at approximately the Rayleigh 
scattering rate which can be written 7r = 
(8/3)E^ki J a2/h. Applying the f-sum rule for 
the dynamic structure factor one can show |19| 
that Rayleigh scattering transfers energy to the 
gas at a rate gj-Etot = ZErNjr. which, sur- 
prisingly, is independent of the interactions be- 
tween the atoms. Comparing this heating rate 
with the energy of the ground state of the gas, 
-Etot ~ H s + fftrap + #kin, where iJ kin is the 
kinetic energy of the atoms, one can estimate 
a heating time via r hcat = E tot / (dE tot / dt) . To 
measure a roton the BEC must survive for longer 



than the roton period r ro ton oc 2ir?i/Ef{ (cf. Fig. 
I2J . For the density and radius stated above, then 
for X = 0.051 one finds Th ca t ~ 8r roton , making 
an experiment challenging but feasible. The sit- 
uation improves for larger, denser, condensates 
since the polarization increases as N 2 , whereas 
7r is a single atom effect. Finally, we note that 
reducing the s-wave scattering length via a Fes- 
hbach resonance allows the laser intensity (and 
hence the Rayleigh scattering) to be reduced by 
an equal factor — see Eq. — and still obtain 
the same effects. 

In conclusion, atom-atom correlations due 
to laser induced dipole-dipole interactions in a 
gaseous condensate can give a roton minimum 
in the Bogoliubov dispersion relation. The cor- 
relations are tunable via parameters such as ra- 
dial width, laser intensity and wavelength. We 
thank M. Boshier, C. Eberlein, J. Steinhauer, 
R. Shiell, and H.T.C. Stoof for illuminating dis- 
cussions, and the Engineering and Physical Sci- 
ences Research Council (EPSRC), the German- 
Israeli Foundation (GIF), and the EU QUACS 
and CQG networks for funding. 
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